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ABSTRACT 

We study the bright end of the luminosity distribution of galaxies in fields with Luminous 
Red Galaxies (LRG) from the Sloan Digital Sky Survey (SDSS). Using 2099 deg 2 of SDSS 
imaging data, we search for luminous (> L*) early-type galaxies within l.Oh" 1 Mpc of a 
volume-limited sample of 12, 608 spectroscopic LRG in the redshift range 0.12 < z < 0.38. 
Most of these objects lie in rich environments, with the LRG being the brightest object within 
1.0/i _1 Mpc. The luminosity gap, M 12 , between the first and second-rank galaxies within 
1.0 hr 1 Mpc is large (~ 0.8 mag), substantially larger than can be explained with an expo- 
nentially decaying luminosity function of galaxies. The brightest member is less luminous 
(by 0.1 - 0.2 mag), and shows a larger gap in LRG selected groups than in cluster-like en- 
vironments. The large luminosity gap shows little evolution with redshift to z — 0.4, ruling 
out the scenario that these LRG selected brightest cluster or group galaxies grow by recent 
cannibalism of cluster members. 

Key words: methods: statistical - galaxies: elliptical and lenticular, cD - galaxies: evolution 
- galaxies: clusters: general 



1 INTRODUCTION 

The bright end of the galaxy luminosity function is still not com- 
pletely characterized. While uncertainties about surface brightness 
completeness is an issue at the faint end iDalcanton 1998), the 
limiting factor at the bright end is often small number statistics. 
Early studies of the bright end of the galaxy luminosity function 
focused on high density regions, often in the richest clusters of 
galaxies. The galaxy number counts drop sharply around the char- 
acteristic luminosity (L*) so fitting functions used to describe 
the luminosity function (expressed in magnitude) have a decay- 
ing bright end - lAbelll ll965l) used a power-law. iHubbld <1936l) 
a Gaussian, while Schechteg ll97q) a double exponential - all 
of which by construction predict a small number at the brightest 
extreme. Hence, the occurrence of even a single galaxy at large 
luminosity (> 10 L*) would seem improbable. However, early 
studies of composite cluster luminosity functions show an ubiqui- 
tous bright end "hump" populat ed by the b rightest members. For 
example, in his seminal paper, Schechter 1 1976) found that the 
Gamma distribution (in luminosity) provides an excellent fit to the 
composite galaxy luminosity function of 13 massive clusters of 
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1 In this paper, we use M* (the characteristic absolute magnitude of the 
knee of the luminosity function) and L„ (the corresponding luminosity) 
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galaxies in the local universe when the brightest galaxy of each 
cluster is excluded from th e fit. This result has bee n reproduced 
by m a ny investigators {Collesj|l98||; |Lu£ge3Ll989|; |yalgtto_et_alJ 
I l997t iLumsden et aljfl997t lTrenthamlll998t iGarilli et alJll999t 
Paolill o et all200ltlYagi et all200AlGoto et all2003l) . 

With the advent of wide-angle redshift surveys, the field 
galaxy luminosity function has been increasing ly well measu red 
iLovedav et alll99llNorberg et all2002llBlanton et all2003 d) . In 
these data, the bright hump seen in the galaxy luminosity func- 
tion derived from rich clusters is either ab sent or not as pronounce. 
When careful correction for EddingtoiJ Jl913l) bias due to pho- 
tometric scatter is included in the likelihood fit to the luminosity 
function, the extreme plunging in number counts predicted from a 
double exponential, is confirmed both locally l Blanton et al. 2003a) 
and up to z ~ 0.4 l Loh 2003). This prompted Schechter 1 2002) 
among others to conclude that the bright end shape of the galaxy lu- 
minosity function is universal- it cuts off exponentially in luminos- 
ity, ma tching the analytical mass function in the lPress & SchechteJ 
1 1974) formalism - lending credence to the mass upper bound cal- 
culations of lRees & Ostrikerl <1977l) . 

In addition to their extreme luminosities, brightest galaxies in 
clusters are often distinguished as having disturbed morphologies 
and extended low surface brightness stellar halos. Ostriker and col- 
laborator s 10Mrikgx&Tre mainell97li : IOstriker & Hausmarlr^77|; 
| Hausman& Ostrikeij|l978l. but see also lRichstonelll975l Il976t 
Iwhitdll976t lBinne\lll97l iGunn & Tinsle^ll977l) ~have proposed 
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that these properties are a consequence of being centrally located 
in high density environments. Indeed, the majority of such objects 
are found in high density environments, and only rarely does the lu- 
minosity function of moderately rich groups (cf. Gell er & Postmanl 
Il983l) exhibit a hump. Hence, accretion via tidal stripping, and 
merger activity via dynamical frinction in rich environments would 
seem to be responsible for their growth. However, careful numeri- 
cal and analytica l modeling of cluster environments (Merritt 1984, 
Il985t iTremainel ["19901 suggest that the luminosities of Brightest 
Cluster Galaxies (BCGs) increase only slightly in a Hubble time, 
as the large velocity dispersions of the galaxy in these rich clus- 
ters makes the proposed cannibalistic scenario inefficient. Canni- 
balism, if it were to operate solely at the present epoch, is insuf- 
ficient to explain the observed magnitude difference Amu (here- 
after the dominance) between the BCG and the second brightest 
galaxy in rich clusters. Recent state of the art numerical simulations 
of massive clusters in a co smological context iGaoetal. 2004; 
lAthanassoula. Garii o&Garcia Gomez 2001; Dubinskvl ll998t see 
also earlier work bv lWesil994l ; lBode et alll9 94) suggest that BCG 
dominance is driven primary by cosmological infall, and it thus de- 
velops as the clusters themselves form. 

In this paper, we us e data from the Sloan Digital Sky Survey 
(SDSS: lYorketalJl2000h to examine the dependence of the lumi- 
nosity difference between the BCG and the second-ranked mem- 
ber on environment and redshift, to test these various evolutionary 
scenarios. Rather than selecting galaxies in a priori chosen dense 
environments, we sample the field for the most luminous red galax- 
ies (i.e., those with the largest stellar masses) at each cosmological 
epoch, using a selection algorithm based on colors and apparent 
magnitudes of a passively evolving old stellar population. 

Specifically, this paper investigates the bright end behavior 
of the luminosity function in regions that host at least one spec- 
trosco pic Luminous Red Galaxy (hereafter LRG; Eisenste in et alJ 
1200 ll) from the SDSS. All high-density regions contain one or 
more LRGs but LRGs also sample lower de nsity environments 
(e.g. lLohetaljE005bl: lEisenstein et"ai] 120051) . including objects 
without companions to the flux limit of the SDSS within ~ 
1.0ft Mpc. This work compleme nts older studies using sam- 
ples selected by galaxy overdensitie s |Tjgmaine& Ri chstoneil977t 
iGeller & Postman 1983; Bhavsar & Barrow 1985) and can test the 
dominance-overdensity hypothesis. It also has a potentially cleaner 
interpretation, as overdensity-selected samples are often subjected 
to bias from projection effects, and by the so-called selection effect 
lScottll957h when comparing samples over a range of redshift. 

The outline of the paper is as follows: In the next section, we 
describe the sample used in our analysis. We describe the statis- 
tical tests used in section [5] In section [4] we outline in detail the 
implementation of our analysis, including a description of the var- 
ious subsamples used. In section|5| we present our results and give 
caveats that need to be taken into account for robust interpretation 
in section [6] Finally, in section [7] we summarize our key conclu- 
sions. Throughout this paper, ft = Ho/100 km s _1 Mpc -1 and 
we adopt the concordance Q m =0.3, Q,a — 0.7 cosmology. 



2 SAMPLE 

2.1 The SDSS Survey 

The SDSS will eventually image a quarter of the Celestial Sphere 
and obtain spectra of approximately 1 million sources. The 5 band 
(ueriz. IFukugita et alJll99rj) imaging is done s imultaneously dur- 
ing photometric conditions iHogg etaljEoOll) using a specially 



designed wide-field camera iGunn et al]|l998l) in drift- scanning 
mode. The imaging data are processed by automated pipelines 
that detect and measure photometric properties iLupton et alJ 
120011 ), and astrometrically and photometrically calibrated the data 
IPieretaljEool; ISmith et aljhool llvezic et aljl2004 . From the 
imaging survey, sources are selected for spectroscopy using a 
640 fiber spectrograph mounted on the same telescope. To date, 
SDSS has had f our major public da ta releases: the Early Data 
Release (EDR: [Stoughton et alJkOo'j) . Data Rel ease One (DR1; 
l Abazaiian et aljbo03h . Data Release Two (DR 2; lAbaz aiian et all 
l2004l) and Data Release Three (DR3: lAbazaiian et alj|2005t) . The 
data used for this paper are drawn from DR1. 



2.2 Luminous Red Galaxy Sample 

The LRG sample I Eisenstein et al. 1200 ll) is a set of intrinsically red 
and luminous (> 3L*) galaxies targeted spectroscopically by SDSS 
to create a large- volume galaxy survey at moderate number density, 
out to z ~ 0.5. The full LRG sample comprises ~ 12% of all 
SDSS galaxy spectroscopic targets IStrauss et alJl2002l) . Here, we 
use those LRGs selected by "Cut I" - galaxies that lie on the linear 
locus of g — r and r — i color space - to give an approximately 
volume-limited sample to rum = 19.2 at z < 0.38. 

The luminosities and colors of giant elliptical galaxies are 
observed to evolve slowly, and aside from the outer stellar enve- 
lope of BCGs and cDs, the stars in these g alaxies are be lieved 
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to have formed a t high redshift, e.g. z > 
Elii^^i2^99jJ^^6^^akmmc^^Mg^^Kauffmann 
Stanford, Eisenhardt & Dickinson 1998; van Dokkum et al. 
Burke, Collins & Mann 2000). The LRG target selection, described 
in detail by |[Eisenstein et alJ 1200 ll) . tunes its selection criteria 
(based on the uniform SDSS photometry) to match the luminos- 
ity and color of a passively evolving old stellar population. The 
K-corrections as a function of redshift are such that the contam- 
ination from intrinsically blue and low-luminosity objects is es- 
sentially negligible for 0.23 < z < 0.38. Morphologically, the 
LRG are bulge-dominated galaxies, and have surface brightness 
distributions and stellar light concentrations similar to those present 
day giant ellipticals and lenticulars. The LRG spectra match the 
spectral energy distribution of an old stellar population, although 
there are unresolved i ssues regarding non-solar abundance ratios 
lEisenstein et ali2 003 ) . 

However, the LRG target selection algorithm only yields a 
consistent population of galaxies for z > 0.23. At lower red- 
shift, the color criteria do not distinguish between less luminous 
and intrinsically bluer galaxies from the desired luminous and red 
ellipticals. However, luminous red ellipticals with z < 0.23 are 
bright enoug h to b e included in the main SDSS galaxy sample 
I Strauss et al. 200 J of which rest frame colors and luminosity to 
select a consistent population. We use an empirically determined 
color-redshift relation, described in Appendix And outlined in 
full in lLoh et all l2005al) . to select z < 0.23 galaxies with the 
same evolution-compensated rest-frame color and luminosities as 
the higher-redshift LRGs. 

FigureQshows the histogram of the comoving density of LRG 
before and after we do the consistent rest-frame selection. Our sam- 
ple has essentially a constant density, and is thus volume-limited 
from z w 0.12 up to roughly z ~ 0.38. The low redshift limit 
is a consequence of the bright-end incompleteness of the SDSS 
spectroscopic survey, as the SDSS spectrographs show excessive 
crosstalk when the flux within the 3" fiber exceeds r ~ 15 (or 
J"pctro ~ 13.5). The high redshift limit comes from the LRG faint 
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Figure 1. The comoving density of LRG (Cut I) as a function of redshift 
before (grey) and after (black) applying the consistent rest-frame color and 
luminosity cut as described in Appendix|A| The new trimmed LRG sample 
is approximately volume-limited to 2 0.38 with a density of ~ 8 X 
W~ 5 h 3 Mpc -3 . 



magnitude limit of rum = 19.2. The comoving density of these 
galaxies is ~ 8 x 10 ~ 5 /i 3 Mpc~ 3 , an order of magnitude larger 
than the abundance of lAbelll 1 19651) Richness Class ^ clusters of 
galaxies, and roughly matches the expected abundance of haloes 
with mass greater than 5 x 10 13 /i" 1 Mq given a concordance 
power spectrum. 

Our primary spectroscopic data are derived from the uniform 
large-s cale structure compilation sample 12 of iBlanton et alJ 
(2004); this covers 2099 deg 2 of sky. A thorough discussio n of 
this sample is given in the appendix of iTeemark et alJ 120041) and 
IBlanton et all COPj) . There are 12, 608 LRG with 0.12 < z < 
0.38 in samplel2 that pass the consistent population rest-frame 
magnitude and color cuts described in Appendix^] We refer to this 
as the trimmed LRG sample. To determine the environments of the 
LRG in this sample, and to measure the dominance of the BCG, 
we use imaging data from the SDSS DR1. We consider all sources 
that are classified as galaxies by the SDSS photometric pipeline 
photo with 7-pctro < 21.0, a level at which the photometry is of 
very high S/N, and at whic h star-gala xy separat i on is particularly 
clean 1 Scranton et al. 2002). r pctro , the Petrosian 1 1976) magnitude 
in the SDSS r-bandpass, is the primar y extended so urce flux mea- 
surement used in this paper; see lStrauss et aljl2002l) for a complete 
discussion. All ma gnitudes and colors are corrected for Galactic 
extinction using the Schleael, Finkbeiner & Davis 1 1998) map and 
assuming Rv =3.1. 



3 STATISTICAL TESTS OF THE LUMINOSITY 
FUNCTION TAIL 

The observ ed narrow luminosity and color distribution of BCGs 
IPostman & LaueJ 19951) . and the possible influence of dynamical 
friction and tidal stripping in their evolution, suggest that their lu- 



minosity function may be distinct from that of the field. This is 
often framed in terms of statistical versus special luminosity func- 
tions: the former argues that BCGs are a mere statistical extreme 
of the luminosity function of cluster ellipticals, while in the latter, 
the BCG forms a different class, having its own unique distribution 
shared among BCGs i n many clust ers. 

Tremaine & Richstone Jl977t) have developed an elegant test 
to distinguish between these two descriptions; this method does not 
require the determination of the galaxy luminosity functions of the 
individual clusters involved. Readers who are interested in the de- 
tailed derivation of the theorems that support the validity of the test 
should refer to their paper. The test hinges on two key assump- 
tions about the nature of the galaxy luminosity distribution in clus- 
ters. The first is lScotd jl957l) s model for the luminosity function, 
which postulates that for a given luminosity distribution, overlap- 
ping magnitude intervals are statistically independent. If ijj is the 
integrated luminosity function 



tp(m) = 



(rn) dm' 



(1) 



then the probability of finding v galaxies in the magnitude interval 

[m a , m b ] is 



p v (m a , m b ) 



[ip(m b ) - i){m a )]" 



exp [tp(m a ) - il)(mbW^ 



The number, v, of galaxies (within a single isolated cluster) is a 
Poisson variate with mean and variance if>(rrib) — ip(m a )- In the 
case of finding a single galaxy in the infinitesimal interval [m,m + 
dm], equation ^2} reduces to 



dP = pi(m,m + dm) = (f>(m)dm 



(3) 



In this model, the probability distribution of the j-th ranked galaxy 
with magnitude rn, pyj (m)dm, is just 2 



py) (m) dm = 



prob {j — 1 galaxies brighter thanm} 
x prob {one galaxy in (m, m + dm)} 
Pj-i (— 00, m) x pi(m, m + dm) 



U - 1) ! 



exp[— ip(m)] c/>(m) dm 



(4) 



Scott's model thus states that galaxy luminosities are drawn from a 
Poisson distribution, with magnitudes m being treated as indepen- 
dent random variables. This is completely embodied in equation 

The second assumption is that the bright end of the luminos- 
ity function drops off as a power-law (exponentially in magnitude), 
which allows one to write down an explicit analytical expression 
for the luminosity distribution of the nth-ranked galaxy. For exam- 
ple, the probability distribution of the luminosity of the first-ranked 
galaxy with a universal bright end differential luminosity function 
4>{m) ~ exp[a(m — mo)] is 



P(i)(m)dm = exp[— tp{m)] 0(m) dm 

~ aexp[a(m-m ) -e Q(m - mo) ]dm. 



(5) 



Here, a parameterizes the steepness of the function, and mo gives 
the normalization (or richn ess) of th e cluster . 

Early studi es by IPeeblesI Jl968l); IPetersonl h970h : 
iGeller & Peebles] j 19761) : ISchechter & Peebles! Il976h that 



2 Numbers in parentheses, e.g. mij\, indicate (j'-th) ranked vari- 
ables. 
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compared observed ranked magnitude distributions with Scott's 
model prediction (equation ^5) and its extension to higher ranks) 
were inconclusive. This is because the catalogs used in these 
studies had insufficient redshift information and reliable cluster 
membership for even the brightest few galaxies. In addition, small 
uncertainties in the relative ^-correction in clusters at various dis- 
tances wash away subtle effects that would otherwise distinguish 
the statistical and special descriptions of the BCG luminosity 
distribution. 



3.1 The Magnitude Difference Between First and 
Second-Ranked Galaxies 



iTremaine & Richstonei ' s key insight was to realize the importance 
of using the magnitude differences for a statistical test. The expo- 
nential cutoff in the differential luminosity function <fi ~ exp(am) 
not only predicts the expected distribution of the first-ranked mag- 
nitudes (equation but also puts a tight constraint on the ex- 
pected distribution of the differences between the ranked mag- 
nitudes. BCGs (the first-ranked members) are observed to have 
a small spread in mag nitude, a (Mi) < 0.3 t Sandagel fl972l 
IPostman & Lauejll995l) . suggesting that on average, the magni- 
tude difference between the first-ranked galaxy, and their respective 
second-ranked galaxies, M12 = M2 — Mi, cannot be too large, if 
both of these gal axies are drawn from th e same exponentially de- 
caying 4>. Indeed. ITremaine & Richstonei showed that the expected 
magnitude difference, (A/12) must at most be of the same order as 
the size of the spread of the first ranked galaxy magnitudes: 



Ti 



a(Mi) 
{M12) 



> 1 



(6) 



Similarly, they showed that the spread of the magnitude difference, 
&(Mi2), must also be of the same order of magnitude as the differ- 
ence (M12): 



o(M 12 ) 
(M12) 



> 0.82 



(7) 



The two inequalities are valid even when there are variations in the 
bright end slope cii of each of the clusters of galaxies used to derive 
the distributions of Mi, M2, M12 etc., so long as the (f>i satisfy the 
two assumptions discussed above. Note that equation J5J holds only 
under the more restrictive assumptions that all a take a single value, 
and that the intrinsic richness of the clusters, mo be the same. 

However, in the limit that all the cti — » (a) < 0, the magni- 
tudes, m, o flScont s model would be an identical and independently 
distributed (i.i.d) random variable. Using extreme- value statistics, 
iBhavsar & Barrowl ll985t) showed that as long as 4> has a gen- 
eral functional form in the domain of attraction of the exponen- 
tial distribution - including the Gaussian, Gamma, exponential and 
double-exponential distributions, then equation |3J is the expected 
distribution for statistical BCG, and the expectation value for our 
two statistics are (Ti) ~ 1.28 and (T 2 ) ~ 1.0. It is remarkable that 
both these values are independent of both a and mo. 

Note that the statistic T 2 depends solely on the differential 
quantity M12, and eliminates uncertainties associated with redshift 
and direction dependent corrections like k+e and extinction, mak- 
ing analyses that combine data from different epochs more robust. 



4 MEASUREMENT OF M u M 2 AND RICHNESS OF LRG 
ENVIRONMENTS 

The trimmed LRG spectroscopic sample we are using extends to 
z — 0.38. However, the main SDSS galaxy sample, probing lower- 
redshift galaxies, has a median redshift of only ~ 0.1, with very 
few galaxies above z = 0.2. We therefore use the photometric 
data from the SDSS to characterize the environments and the lu- 
minosity distribution of galaxies around each LRG, and to calcu- 
late statistics such as M12. We do so by using color cuts to iso- 
late early-type galaxies at the same redshift as the LRG, using the 
red sequence seen in cluster color-magnitude diagrams. Within an 
angular extent equivalent to 1.0 hT x Mpc at the redshift of each 
spectroscopic LRG, we searched the SDSS photometric data for 
second- and third-ranked galaxies with g — r and r — i colors typ- 
ical of early-type galaxies at that redshift. The regularity of the ob- 
served color-magnitude relation of early-type g alaxies seen out to 
z ~ 0.5 and beyond, e.g. Blakesl ee et alJ 120031) . enables us to em- 
pirically determine a color-redshift locus for these galaxies. This is 
shown in Figure IbTI with an estimated dispersion at each redshift 
given by the solid ellipses. Details of how this locus is determined 
is described in AppendixlBl Note that the locus is not simply drawn 
from the colors of LRG themselves since "Cut I" LRGs are red bi- 
ased - lower luminosity galaxies need to have a redder color to pass 
the selection cut (Eisenstein et al. 2001). The ellipses include the 
observed color errors - almost entirely due to photometric scatter 3 
- of galaxies on the red sequence with apparent magnitude brighter 
than ~ M* (z) at each redshift. Hence, if we include only galax- 
ies from the fields around each LRG with colors within the n a 
ellipse of the LRG redshift, then on average, we will be complete 
with respect to red sequence galaxies in the LRG fields at the n a 
level. Of course, for larger n, the background contamination goes 
up. We find that the contamination remains manageable, while still 
including the majority of galaxies in the red sequence, if we include 
galaxies that are within 2.5 a of the color locus. 

The galaxies are then ranked by apparent Petrosian r magni- 
tude. The first-ranked galaxy in each field is the LRG 95-99% of 
the time, with the lower values holding in higher redshift fields. 
In some of the exceptions, the first-ranked galaxy was also flagged 
as an LRG candidate, but a spectrum was not obtained because of 
the restric tion that two spectros copic fibers not be placed closer 
than 55" iBlan ton et alj|200 3a). The majority of the remainder 
turned out to be bluer galaxies presumably in the foreground, se- 
lected because of the generosity of the 2.5 a color cut. We keep 
this small number of non-LRG first-ranked galaxies in our analysis 
so as not to bias our result against the occasional true bluer first- 
ranked galaxy. However, excluding them from our sample does not 
significantly change our results. 

Our first goal is to derive the intrinsic joint distribution of Mi 
and M12, 4>(Mi, M12). For this, we need to estimate and subtract 
the background distribution <^ back (Mi, M12). For each LRG field, 
we select ten random posi tions within the same SDSS imaging 
stripe iYork et aT]|2000l) . Just as in the LRG fields, we draw 
a 1.0 h Mpc circle, apply the same 2.5 a red sequence color 
cut, and rank the galaxies by brightness. The <j> are binned with 
AMi — 0.1 mag and AA/12 = 0.2 mag. Figures |2| and [3] show 
the composite <j> for 271 and 1732 LRG fields at 0.14 < z < 0.16, 
and 0.34 < z < 0.36, respectively. The top left panel of each 



3 The intrinsic scatter of the color-magnitude relati on is expected to be less 
than 0.05 magni tude even at z = 0.4 lEUis et alj|l997tlCool et al]|2005t 
lLoh et all2005al) . 
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Figure 2. These figures illustrate the background subtraction in Afi and Afi2 space, done in a narrow redshift interval of 0.14 < z < 0.16. The top left 
panel gives the observed composite <f>(Mi, M12) for 271 LRG fields; top right is the background estimated from 3210 random positions, but scaled by area to 
match the observed <j>(Mi , M12). The bottom left panel is the corrected distribution using equation |S}, weighting by / = 228/271, the ratio of the number 
of uncontaminated (effective) fields to the total observed fields. The bottom right gives the marginalized 1 st ranked distribution. The blue histogram is the 
observed (j>(Mi) distribution, and the black histogram is the effective distribution obtained by subtracting the (dashed) background. The solid vertical line 
gives the mean (Mi). The magnitude of a M* — 1.5 early-type galaxy at the median redshift of 2 ~ 0.15 is indicated by the dashed vertical line in all panels. 
Density contours and histograms are normalized to unity. 



figure is the observed distribution while the bottom left gives the 
background. The dashed vertical line is the observed magnitude of 
a M* — 1.5 galaxy at the respective redshift, while the slanted solid 
line gives the upper bound to the Mi 2 we could observe given the 
SDSS photometric catalog limit of r pc t ro = 21.0. Betwee n 0.5% 
and 1% of all LRG in our sample are isolated fsee lLoh et all2005bl 
for further discussion), showing no companion within the color and 
magnitude space we searched. These objects may well be fossil 
groups IPonman et alll994l) . showing substantial numbers of low- 



luminosity companions below our photometric limits. We leave the 
study of these objects to a future paper. 

We estimate the intrinsic distribution in Mi and M12 (given 
in the top right panels) using a weighted subtraction given by 

$ (Mi , Mi 2 ) = /x [0 obs (Mi , Mia) - a< back (Mi , M 12 )] 

+(1-/) x0 obs (Mi,Mi 3 ) (8) 

where / = ( jV bs — Nb&ck)/N bs is the fraction of fields in which 
M12 measurements are not affected by contamination, and N bs 
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Figure 3. As in Fig. [2] for 0.34 < z < 0.36. The solid diagonal line in the 2D plots gives the upper bound of M\i for a given M\ set by the flux limit 
(fpctro < 21.0) of the SDSS imaging catalog. 



and iVback are the number of LRG and background fields in the 
Mi,Mi2 bin in question. If the background contamination (1 — 
/) is large, there is a substantial probability that the galaxy we've 
labelled the second-ranked is actually a background (or foreground) 
object, thus what we have labelled the third-ranked galaxy is in fact 
M2. This is the origin of the second term in equation {8}. 

We first perform our analysis on the full volume-limited 
trimmed sample of LRG in a series of redshift intervals of width 
Az = 0.02 from z = 0.12 to z = 0.38. We define the richness (or 
the associated galaxy overdensity) of each LRG field by counting 
the number of galaxies on the red sequence withi n 1.0 fa" 1 Mpc of 
the LRG from M* - 2.5 down to M* lLoh et all2005bT) . We then 
split the sample into four, in terms of quartiles based on the rich- 
ness of these LRG fields; the richness ranking is done within each 



redshift slice. We term the relative rank of each LRG in its red- 
shift slice the .R-level. The upper quartile sample, with a median 
of 8 early-type galaxies with magnitude M* and brighter within a 
1.0 Mpc radius, are th e equivalent o f a present day moderately 
rich cluster of galaxies (e.g. Abell ( 1965) Richness Class > 1). Ob- 
jects in the 25 th — 50 th percentile, on the other hand, with 2 to 3 
galaxies on the red sequence, are similar to present day groups of 
galaxies. 

Figure[4]is the redshift histogram of LRG used for our anal- 
ysis. The blue shading refers to the upper quartile, while the red 
shading indicates the 25 th — 50 th percentile. Note that the fraction 
of LRG fields in the upper quartile (say) of a given redshift sub- 
sample can be greater or smaller than 25% of the subsample total 
because LRGs with the same richness share the same rank. 
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Figure 4. Redshift histogram of the trimmed LRG sample. The blue (red) 
shading indicate the number of LRG fields with richness in the upper 
75 th (25 th - 50 th ) percentile. 



5 TREMAINE-RICHSTONE STATISTICS AS 
FUNCTIONS OF REDSHIFT AND RICHNESS 

Figures[2]and[3]show that as one goes to higher redshift, LRG clus- 
ters tend to have smaller Mi 2. For example, at 2 ~ 0.15, there are 
quite a few LRG clusters with M12 ~ 2, well above the noise level 
(given by the bottom left panels), while there are no clusters with 
such a large gap at z w 0.35. 

This systematic smaller gap at higher redshift is also seen 
when the sample is split into group-like (25 th — 50 th ) and cluster- 
like (upper quartile) LRG fields, as shown in the top panels of Fig- 
ures|5|and|6| The bottom panels of each figure show histograms of 
the marginal gap distribution, 0(Mi2), explicitly showing the noise 
level, back (Mi2) (dashed) determined from 10 random positions 
for each LRG, the observed <p° bs (Mi 2 ) (blue) and (t>° bs {M 13 ) 
(red), and the final distribution 0(Mi2) (solid) used for our analy- 
ses as estimated from equation {8j. The vertical dashed line is the 
mean (M12) for each panel, estimated from the first moment of 
(/>(Mi2). It shows a slight decrease with redshift for both clusters 
and groups. 

The iBautz & Moreanl ll970l) effect, whereby clusters with 
brighter BCGs have larger gaps, is seen in both rich and poor sys- 
tems at all redshifts. One measure of the strength of this effect is 
the tilt of the number density ellipsoid of the joint distribution, 
cj>(Mi, M12). From Figures [5] and |6| it appears that this effect is 
stronger in poorer systems, but changes little with redshift. 

LRG fields with group-like environments have a larger gap 
than do cluster-like environments at all redshifts. This is seen in 
both the joint distribution, <f>(Mi, M12), and marginal distribution, 
0(Mi2) (Figures [5] and This must in part be statistical: Pois- 
son sampling from a common luminosity function naturally gives a 
smaller gap for a larger total number of galaxies. We will see below 
from the Tremaine-Richstone statistics that this is not the whole 
story; Figure [8] shows that the bounds of equations |S| and Q are 
violated. If the universal galaxy luminosity function estimated over 



large scale (cf. iBlanton et alJ|20 03al) is correct, i.e. the bright end 
cuts off exponentially like"a ]schechteJ function, then the Poisson 
assumption must not hold. 

FigureQsummarizes the evolution of the mean and dispersion 
of Mi and M12 from the various sub-samples. The top right panel 
shows that poorer systems have first-rank members that are brighter 
by 0.1-0.2 mag in the mean, while the top left panel shows that 
these same systems have a larger mean (M12) at all redshifts. Our 
sample is volume-limited independent of richness, and the second 
and third ranked galaxies are also complete for the range of A/12 
and M13 we are probing, given the depth of the SDSS photometric 
data. The dependence of the gap size on richness is not an artifact 
of our richness estimates, as the mean (and median) number of M* 
and brighter galaxies is roughly constant with redshift, after back- 
ground subtraction (cf. lLoh et all2 005b). 



Our results for the iTremaine & Richstond statistics Ti and 
T2 are presented in Figure [8] The estimated Ti and T2 are well 
below the lower bounds for the statistical BCG hypothesis (indi- 
cated by the solid horizontal line) set by equation {6j and Q for 
both poor and rich systems at all redshifts. They are even fur- 
ther away from the statistical i.i .d. prediction (uppe r dashed line) 
from extreme-value statistics of Bha vsar & Barrowl ^985). If we 
assume that the bright end of the galaxy luminosity function in 
clusters cuts off exponentially like the universal galaxy luminos- 
ity function, then the statistical hypothesis for BCG is rejected with 
high confidence, even if we allow the steepness of this cut-off to 
vary from cluster to cluster. Of course, this result has meaning 
only if our initial assumption about the galaxy luminosity func- 
tion embodied in lScottl 's model - of independent luminosity sam- 
pling - is valid. The assumption that the luminosities of galaxies 
are independent must break down at some level, since we know 
that galaxies influence with each other via tidal interactions and 
mergers through gravity to very large distances. On the scale of 
1.0 Mpc, the independence assumption is unlikely to be true 
especially for galaxies as massive as LRGs (~ 5 - 8 x 1O 12 M ). 
Hence the fact that Ti and T2 fall well below the boundaries of 
equations {5J and Q does not by itself imply that BCGs follow 
an assembly path different from that of other cluster galaxies. It 
might merely signal a breakdown in our assumption about the 

galaxy luminosity function. Indeed, bright elliptical galaxies, in- 

^^^^^^^^^^^^^^^^^^^^^^^ 

eluding BC Gs 1 Oeger le & Hoessellll99lt Postman & LaueJll995t 

iGraham et aljll998l) for m a low dimensiona l sequence (the Fun- 
dame ntal Plane relation, iDressler et alJll987t iDiorgovski & Davisl 
Il987t) . implying some kind of universal assembly history. 



Ti and T2 are the inverse of the mean luminosity ratio of the 
two brightest galaxies, measured in units of the two characteris- 
tic luminosity spreads, cr(Mi) and c(Mt2). When considered as 
a function of redshift, they quantify the normalized relative bright- 
ening of the first-ranked galaxy, as compared with lesser members 
of the same cluster. These normalizations are with respect to the 
dispersions of Mi and M12. Ti(z) and 22(2) are consistent with 
being flat as a function of redshift, suggesting that there is no ma- 
jor process that influences the creation of the observed luminosity 
gap over the redshift range we probe. This can be compared with 
the raw gap function, (Mi2)(;z), of the top left panel in Figure^ 
where a decrease with redshift is seen, suggesting mild brightening 
of the first-ranked galaxy compared to the second-ranked galaxy 
with time. Nonetheless, the corresponding spread in magnitude, no- 
tably a (A/12), also decreases with redshift, in such a way that T2 
remains flat. 
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Figure 5. On the top left (right) panel is <f>(Mi , M12) for LRG fields with richness .R-level in the 25 th to 50 th (75 th and above) percentile. Higher percentile 
corresponds to richer systems. The bottom panels are the corresponding marginalized 4 } {Mi2)- The black histogram is the final distribution estimated from 
equation jS} - a weighted linear combination of the observed Mi 2 (blue) and (M13) (red) distribution, and the background distribution (dashed). For LRG 
fields with redshifts 0.14 < z < 0.16. 



6 DISCUSSION 

6.1 Interpretive Caveats 

Our richness measure says little about the gravitational potential 
of the LRG fields, since we do not measure the mass overden- 
sity directly. Our proxy for overdensity, the .R-level yardstick, is 
tuned to robustly characterize the environment of rich systems with 
high signal-to-noise ratio. Our groups (i.e., the 25* — 50 th per- 
centile sample) by design must host a luminous (< A/* — 1 or 
brighter) LRG, hence it is a highly biased sub-class of galaxy 
groups. The comoving number density of our groups is a mere 



2 x 10" 5 /i 3 Mpc- 3 , an order of magnitu de below the exp ected 
group density from the analytical IPress & SchechteJ J 1974) esti- 
mates. Hence, we do not sample a substantial fraction of group-like 
environments in the universe. Even if the mass -to-iiMevel scaling 
relation is tight and unbiased <Loh et all2005cK our richness mea- 
surements can give biased mean mass estimates in the present of 
moderate intrinsic scatter and measurement uncertainties. Further, 
it could be that our groups are equally massive, with a huge dwarf 
galaxy population compared to their richer counterparts, as ob- 
served loca lly in a subset of X-ra y bright Hickson compact groups 
jHunsberger. Charlton & Zaritskvll99ot . 
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Figure 6. As in Fig.|5] for redshifts 0.34 < z < 0.36. 



These groups are h i gher r edshift analogs of the famous 
iMorgan. Kavser & White! 1 19751) and lAlbert White & Morgan! 
1 1977) groups that host luminous (non-c) D galaxies. However, 
we are only sampling a subset of special groups, which may re- 
solve t he apparent contradiction between our results and those of 
iGeller & Postmarl Jl983l) - they found that the brightest members 
of groups are less dominant than those of clusters and are consis- 
tent with them simply representing the luminous tail of the lumi- 
nosity function. It would be of interest to know what fraction of all 
groups (observationally defined) host an LRG. We will show in a 
subsequent paper <Loh et all2 005c> that ~ 95% of the most mas- 
sive clusters selected by both the optical matched-filter techniques 
<Kim etall2 002a) and X-rays have LRGs in them, but the fraction 
for groups is expected to drop substantially. Note that the num- 



ber density of groups from optical redshift surveys have a higher 
number density than LRGs, in part because many of these groups 
are not actually virialized and display no extended X-ray emission. 
Further, there are groups with extended X-ray emission that have 
bright elliptical members that are just shy of the M* — 1 luminos- 
ity threshold of the LRG <Nipoti et all2003h . 

Another caveat of our analysis is the 1.0 h~ x Mpc scale used 
to search for second and third ranked galaxies. If we choose a 
smaller scale, say 0.25/i _1 Mpc, we find a larger mean gap for 
both richer and poorer system, by ~ 0.5 mag. However, the dif- 
ference in the gap between rich and poor systems is smaller. For 
example, at the current l.Ofo -1 Mpc scale, the average difference 
of the estimated (M12) between the rich and poor systems is 0.4 
mag (see the top right panel of Figure^, while at 0.25/i _1 Mpc 
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Figure 7. (Left Panels) Average 1 st ranked r pe t ro magnitude (top) and its corresponding 1 a dispersion (bottom) as a function of redshift. (Right Panels) 
(M12) and its corresponding 1 a dispersion as a function of redshift. The blue (red) is for fields with richness in the upper 75 th (25 th - 50 th ) percentile. 
EiTors bars are lcr from bootstrap. The dot-dashed line is at z = 0.36, the redshift beyond which the sample is no longer volume-limited. 



it is ~ 0.15 mag. This difference is again in part Poissonian, but 
is also influenced by how bright galaxies populate a cluster. The 
smaller scale analysis of 0.25/i -1 Mpc would, for example, split 
the Coma cluster into two LRG fields, each with a large luminosity 
gap, while the present l.O/i -1 Mpc analysis would treat Coma as a 
single system with a small gap. We could try to infer a dynamically 
motivated scale like a virial radius or V2oo - the radius at which the 
galaxy density drops below 200 times the mean density - but these 
measurements are noisy and ill defined for poor systems. 



6.2 Comparison with Numerical Studies 

We compare o ur findings with the numerical experiment of 
iDubinskvl ll998t) . who simulated the formation of groups and poor 
clusters of galaxies in a cosmological setting (albeit based on 
the now out-of-favor Q m — l,as = 0.7 CDM model). Our 
25 th — 50 th percentile i?-level sample matches his systems. His 
simulation specifically addresses the growth of the brightest galaxy 
and its creation in the context of cluster collapse in a hierarchi- 
cal structure formation paradigm. His experiment begins at 2 = 2, 
and by z ~ 0.8 the major galactic unit that is later identified as 
the BCG is already in place. This BCG is a result of mergers of 
the brightest few galaxies of roughly equal mass during the triaxial 
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Figure 8. The Tremaine & Ri chstond statistics T\ and T2 as a function of redshift computed in redshift slices Az = 0.02. The solid horizontal line 

gives bounds below which brightest members would cease to be the statistical extreme of the luminosity function of lesser members in the cluster. The top 
dashed line is the prediction from extreme- value statistics. The blue (red) curve is for fields with richness in the upper 75 th (25 th — 50 th ) percentile. Errors 
bars are la from bootstrap. The dot-dashed line is at z = 0.36, the redshift beyond which the sample is no longer volume-limited. 



collapse of the cluster. By z ~ 0.4, all major mergers between the 
BCG and other galactic units are complete. The large mass contrast 
between the BCG and the second-ranked galaxy is now in place. 
This would suggest that M12 will shrink at z > 0.4, unfortunately 
beyond the redshift which our sample probes. Qualitatively, our re- 
sults from the group-like LRG fields are in agreement with his nu- 
merical experiment. The large and unchanging gap seen in these 
poorer systems points to an earlier origin of the growth of the dom- 
inance of the BCG, of which the picture painted by his experiment 
is one plausible description. However, cluster infall and growth is 
a strong function of Q. m and erg (in particular, models with smaller 
Sl m predict infall ending at much higher redshift). Thus there is an 
urgent need for simulations of BCG growth using the concordance 
cosmology. 



6.3 Are First-Ranked LRGs drawn from the Universal 
Luminosity Function? 

BCG show a large luminosity gap; moreover, the dominance of 
the BCG is observed to be independent of redshift. Thus sug- 
gests that whatever mechanism created this dominance must have 
taken place at z > 0.4. Scenarios in which recent c annibalistic 
activity fueled the BCG growth (e.g. lHausman & Ostrike3ll978h 
are unlikely. Poorer systems have more dominant BCGs, suggest- 
ing that the process that makes these galaxies unusually luminous 
is not related to the presence of BCGs in crowded rich cluster 
cores. U sing 211 c l usters o f galaxies selected from SDSS pho- 
tometry, iKim et alj J2002tl) found that BCG with M12 > 0.5 
mag showed strong alignment with their host clusters, while those 
with weaker dominance show ed no tendency to align. (See also 
iFuller. West & Bridged 1 19991) for similar results for MKW/AWM 
group s.) Taken to gether, these observations support the scenario of 
Dubinskv 1 1998) - essentially a generic feature of the hierarchi- 



cal structure formation paradigm with Gaussian initial conditions - 
whereby major mergers during galactic infall along filaments create 
the BCG at moderate (2 ~ 0.5 — 1.0) redshifts. Statistical sam- 
pling may explain why poorer systems (which host an LRG) have 
more dominant brightest members. Alternatively, these systems are 
preferentially located within older structures, and are dynamically 
more mature. The BCGs at the confluence of these structures are 
further along the evolution sequence, making them more dominant. 
Nevertheless, it suggests that the properties of BCGs are intimately 
related to the assembly history of the local matter density. These 
findings, however, run contra ry to the stan dard candle nat ure of 
the BCG metric lumino sities JSandagdll972l iGunn & Okdll975t 
IPostman & Lauejll995l) . and argue against a BCG special popu- 
lation independent of local conditions. Properties of BCG clearly 
depend on environment and yet they are standard candles. We do 
not have a full understanding of this dichotomy. 



7 CONCLUSIONS 

We have shown that the luminosity gap, M12 between the first- 
ranked galaxy and the second-ranked in fields containing at least 
one LRG is large, and this is inconsistent with these brightest 
members being the statistical extre me of the local galaxy lumi- 
nosity function within Scotg ll957l) 's model. This dominance of 
the brightest members changes little over the redshift range of 
0.12 < z < 0.38, and is more prominent in the poorer group- 
like LRG fields than the richer cluster-like fields. Furthermore, 
the brightest member in cluster-like LRG fields are systematically 
brighter (0.1-0.2 mag) in the mean than in group LRG fields. 
When w e co mbine these findings with the c onclusions of Ki m et all 
<2002bl) and lFuller. West & Bridget (l999) that dominant brightest 
members are preferentially aligned with their host clusters, our re- 
sults support the assembly of BCGs during cluster collapse within 
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the framework of hierarchical structure formation, and put strin- 
gent constraints on the growth of BCGs through recent accretion 
of lesser members. In future papers, we will present a statistical 
analysis of the o ptical env ironments of this volume-limited sample 
of LRG jLoh et alfeOO Sb), and their relationship with intermediate 
redshift clusters selected by X-ray emission lLoh et all2005cl) . 
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APPENDIX A: A CONSISTENT POPULATION OF LRG 

As discussed earlier in the main text and in Eise nstein et alj 1200 ll) . 
the LRG target selection will only give a consistent population for 
2 > 0.23. In particular, for r < 17.5, the selection admits many 
intrinsically fainter galaxies. Furthermore, bright galaxies are at 
the exponential tail of the luminosity function, and small devia- 
tions from the nominal LRG luminosity threshold affect the derived 
number counts substantially. Since the flux-limited MAIN galaxy 
survey is spectroscopically complete to r ~ 17.8 or M* — 1 at 
2 ~ 0.25, we spectroscopically trimmed the low redshift LRG to 
match the rest-frame color and luminosity of their higher redshift 
counterparts. We use the empirically derived color-redshift relation 
of M* — 1 red sequence galaxies from Appendix IbI to extract a 
co-evolving population of similar color and luminosity. 

Figure IA1I shows plots of the color-magnitude relation for 
galaxies from the SDSS MAIN galaxy survey. The color used is 
clrg = 7/9g + 5/3r — 4i, which is approximately parallel to the 
color locus of LRGs; see Figure lBll The figure shows three narrow 
redshift slices (Az = 0.01) at z = 0.10, 0.15, and 0.20. The red 
sequence is clearly visible. The axes for the first three panels (both 
top and bottom left) are labeled in observed quantities, but shifted 
relative to the red sequence for direct comparison between the red- 
shift subsamples. The fourth panel (bottom right) is a combined 
plot of the three, where the axis is now labeled in rest-frame quanti- 
ties. We employ a single model k+e correction of a passively evolv- 
ing early-type M* galaxy from PEGASE to realign the r pct ro-axis 



for rest-frame comparison. Galaxies that pass the LRG target selec- 
tion are indicated by violet points. The steep oblique selection edge 
is the color-luminosity bound that imposed a red bias to LRGs. No- 
tice that the selection is more permissive at z = 0.1, but becomes 
gradually more stringent at z ~ 0.2 4 . At z — 0.1 (top left panel), a 
~ M* — 0.5 galaxy with the exact red-sequence color (i.e. one that 
lies on the solid line) would pass the LRG target selection (and be 
colored violet), but at z — 0.2 (bottom left), it would have to have 
an absolute magnitude ~ M* — 1 to make the cut. The red box 
is the proposed rest-frame color-luminosity cut. When applied uni- 
formly across cosmological epoch, we obtained a trimmed sample 
of consistent LRG population with a roughly constant comoving 
density (cf. Figure^. The slight increase in density at higher red- 
shift is not real but a consequence of Eddington bias from larger 
photometric uncertainties at high redshift. 



APPENDIX B: EMPIRICAL COLOR-REDSHIFT 
RELATION 

Since galaxies are clustered, their neighbors tend to be at the same 
redshift. Hence, by using a set of reference galaxies (LRG or QSO) 
with spectroscopic redshifts, one can infer the average unbiased 
color (and other photometric attributes) of a typical galaxy at those 
redshifts from the properties of galaxies that clustered around these 
reference galaxies, even though the selection of the original refer- 
ence galaxies maybe biased. Specifically, for spectra obtained with 
the SDSS, while both LRGs and QSOs are color selected, and thus 
by design have a biased mean color, the more numerous galaxies 
from the deeper imaging data that cluster around these spectro- 
scopic objects are not. By stacking fields of LRGs or QSOs with 
similar redshifts to make a composite galaxy field at those redshifts, 
background and foreground galaxies seen in projection can be re- 
moved statistically as they are not expected to correlate with the 
reference spectra. 

We have calibrated the color-redshift relation and the observed 
scatter, in two color dimensions (g — r and r — i) for galaxies with 
a fixed absolute magnitude range (M* — 2.5 < r pc t r o < A/*). An 
evolving model of an old stellar population elliptical galaxies from 
PEGASE (Fioc & Rocca-Volmerange 1997) was used to compute 
the r-band magnitude for a M* galaxy as a function of redshift. 
Two luminosity ranges were computed: (1) using M* — 2.5 < 
r pct ro < M* - 0.8 for ~ M * - 1 LRG-like galaxies, (2) M* - 
1.5 < r pc tro < M* + 0.5 for ~ M* red sequence galaxies. We 
employ an optimal weighting using a two dimensional elliptical 
Gaussian to locate the peak and estimate the scatter of the locus. 
Details on how this is done and measurements of intrinsic scatter 
of the red sequence to constrain the spread of formation time scale 
of old stars as a function of redshift, as well as comparison between 
colors of galaxies in the neighborhood of LRGs a nd QSOs will be 
discussed in a subsequent paper <Lohetall20 05a). 

Figure iBTl is the g — r and r — i color-color diagram of red- 
sequence galaxies with r pc tro ~ M* , The solid ellipses are the 

4 As an interesting aside, notice that even though the LRG target 
selection is permissive at z « 0.10, it still primarily selects galax- 
ies on the red-sequence, as contamination from blue star-forming 
galaxies is minimal. Thus one could in principle cleanly extend 
the LRG selection of red galaxies to ~ 0.8 mag lower luminos- 
ity, inc reasing the LRG comoving number density by a factor ~ 3 
( Padmanab han et all2004l) 



14 



Y.-S. Loh and M. A. Strauss 




Figure Al. These panels show the observed color magnitude diagrams of LRG and MAIN sample galaxies in three redshift ranges. The axes for the two top 
and the bottom left panels are in observed quantities. A model k+e correction for an old stellar population sequence is used to shift the x-axes to the same 
luminosity range. M* at the respective redshifts of (0.1, 0.15, 0.2 (Az = 0.01)) of the panels are aligned, and is indicated by the dashed vertical line. The 
location of the red sequence - the sloping horizontal line - for a fixed absolute magnitude is used to shift the y-axis to the same rest-frame color. The diagrams 
are made using spectroscopic galaxies brighter than r pctro ~ 17.8 from the MAIN sample. Hence, they are unbiased with respect to color. The violet points 
are these MAIN galaxies that also pass the LRG Cut I selection criteria. The red box is the rest-frame color-magnitude cut used here for the consistent trimmed 
LRG sample. The last panel (bottom right) plots all galaxies from the other three panels, now shown in rest-frame quantities. 



lcr contours while the gray ellipses are 3o\ The locus of early type 
galaxy follows an approximate linear sequence with redshift up to 
the up-turn point at around z 0.37, m otivating the "Cut I" of the 
LRG target selection (Eise nstein et all2 0011. There is a slight kink 
at z ~ 0.26, which corresponds to the M star locus crossing with 
this early-type locus. The analysis was done using ~ 30, 000 LRG 
and QSO in the redshift range of 0.08 < z < 0.44 and 2340 deg 2 
of DR1 imaging data processed with photo 5.3 photometry. 
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Emripical color— redshift relation 
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Figure Bl. The empirically determined locus of red sequence M* galaxies 
from z = 0.09 to z = 0.43 in increments of Az = 0.02. The ellipses 
in black (grey) are the 1 (3) a dispersion of colors at each redshift. This 
sequence of ellipses serves as a pseudo photometric redshift filter to reduce 
contamination from background galaxies. 



